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Abstract 

We introduce a new class of frames with strong symmetry properties called geometrically 
uniform frames (GU), that are defined over an abelian group of unitary matrices and are gen- 
erated by a single generating vector. The notion of GU frames is then extended to compound 
GU (GGU) frames which are generated by an abelian group of unitary matrices using multiple 
generating vectors. 

The dual frame vectors and canonical tight frame vectors associated with GU frames are 
shown to be GU and therefore generated by a single generating vector, which can be computed 
very efficiently using a Fourier transform defined over the generating group of the frame. Simi- 
larly, the dual frame vectors and canonical tight frame vectors associated with CGU frames are 
shown to be CGU. 

The impact of removing single or multiple elements from a GU frame is considered. A 
systematic method for constructing optimal GU frames from a given set of frame vectors that 
are not GU is also developed. Finally, the Euclidean distance properties of GU frames are 
discussed and conditions are derived on the abelian group of unitary matrices to yield GU 
frames with strictly positive distance spectrum irrespective of the generating vector. 
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1 Introduction 

Frames are generalizations of bases which lead to redundant signal expansions |Q, §]• A (finite) 
frame for a Hilbert space W is a set of not necessarily linearly independent vectors that spans 7i. 
Since the frame vectors can be linearly dependent, the conditions on frame vectors are usually not 
as stringent as the conditions on bases, allowing for increased flexibility in their design |^ 

Frames were first introduced by Duffin and Schaeffer in the context of nonharmonic Fourier 
series, and play an important role in the theory of nonuniform sampling ||l|, ^, |[ and wavelet theory 
[P, 1^. Recently, frames have been used to analyze and design oversampled filter banks [0, |8|, ^ and 



error correction codes [10]. Frames have also been applied to the development of modern uniform 



and nonuniform sampling techniques |11], to various detection problems [12, 13 [, and to multiple 
description source coding [0]. 

Two important classes of highly structured frames are Gabor (Weyl-Heisenberg (WH)) frames 
[^, [l^ and wavelet frames [^, ^ Both classes of frames are generated by a single generating 
function. WH frames are obtained by translations and modulations of the generating function 
(referred to as the window function), and wavelet frames are obtained by shifts and dilations of the 
generating function (referred to as the mother wavelet). In Section ^ of this paper, we introduce 
a new class of frames which we refer to as geometrically uniform (GU) frames, that like WH and 
wavelet frames are generated from a single generating vector. These frames are defined by an 
abelian group Q of unitary matrices, referred to as the generating group of the frame. We note 
that WH frames and wavelet frames are in general not GU since the underlying group of matrices 
is in general not abelian. GU frames are based on the notion of GU vector sets first introduced by 



Slepian [18[ and later extended by Forney [|^, which are known to have strong symmetry properties 



that may be desirable in various applications such as channel coding [^, 2C, 21 [. 

The notion of GU frames is then extended to frames that are generated by an abelian group 
Q of unitary matrices using multiple generating vectors. Such frames are not necessarily GU, but 



consist of subsets of GU vector sets that are each generated by Q. We refer to this class of frames 
as compound GU (CGU) frames, and develop their properties in Section |6|. CGU frames are a 
generalization of filter bank frames introduced in [@, |9|. An interesting class of frames results 
when the set of generating vectors is itself GU, generated by an abelian group Q. (Note that this 
class of frames will in general not be GU). As we show, these frames are a generalization of WH 
frames in which Q is the group of translations and Q is the group of modulations. 

Given a frame for TC, any signal in Ti. can be represented as a linear combination of the frame 
vectors. However, if the frame vectors are linearly dependent, then the coefficients in this expansion 
are not unique. A popular choice of coefficients are the inner products of the signal with a set of 
analysis frame vectors called the dual frame vectors [^] . This choice of coefficients has the property 
that among all possible coefficients it has the minimal Z2-iiorm |0, 22 1. 

In Section we show that the dual frame vectors associated with a GU frame are also GU, and 
therefore generated by a single generating vector. Furthermore, we demonstrate that the generating 
vector can be computed very efficiently using a Fourier transform defined over the generating group 
Q of the frame. Similarly, in Section ^ we show that the dual frame vectors associated with a CGU 
frame are also CGU. When the generating vectors of the CGU frame are GU and generated by a 
group G that commutes up to a phase factor with the group Q, the dual frame is generated by a 
single generating vector, a result well known for WH frames. 

An important topic in frame theory is the behavior of a frame when elements of the frame are 
removed. In Section |^, we show that the frame bounds of the frame resulting from removing a single 
vector of a GU frame are the same regardless of the particular vector removed. In this sense GU 
frames exhibit an interesting robustness property which is of particular importance in applications 



such as multiple description source coding |14]. We also consider the behavior of a GU frame when 
groups of frame elements are removed. 

In the special case of a tight frame the dual frame vectors are proportional to the original 
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frame vectors so that the reconstruction formula is particularly simple. In many applications it is 
therefore desirable to construct a tight frame from an arbitrary set of frame vectors. A popular 
tight frame construction is the so-called canonical tight frame []T7|, |7|, |2^, 24, 25, 26 1, first proposed 
in the context of wavelets in [27|. The canonical tight frame is relatively simple to construct, it is 



optimal in a least-squares sense 26, 29 1, it can be determined directly from the given vectors, 
and plays an important role in wavelet theory pO| , 31, 32 1. Like the dual frame vectors, we show 
that the canonical tight frame vectors associated with a GU frame are GU, and the canonical tight 
frame vectors associated with a CGU frame are CGU. When the generating vectors of the CGU 
frame are GU and generated by a group Q that commutes up to a phase factor with Q, the canonical 
tight frame vectors can be obtained by a single generating vector, generalizing a result well-known 
in WH frame theory. 

Since GU frames have nice symmetry properties, it may be desirable to construct such a frame 
from a given set of frame vectors. The problem of frame design has received relatively little 
attention in the literature. Systematic methods for constructing optimal tight frames have been 
considered |25, 2(:, Methods for generating frames starting from a given frame are described 
in iQ]. In Section ^, we systematically construct optimal GU frames from a given set of vectors, 
that are closest in a least-squares sense to the original frame vectors. We consider three different 
constraints on the GU frame vectors. First, we treat the case in which the inner products of the 
frame vectors are known. The optimizing frame is referred to as the scaled-constrained least-squares 
GU frame (SC-LSGUF). Next, we consider the case where the inner products are known up to a 
scale factor. The optimizing frame in this case is referred to as the constrained least-squares GU 
frame (G-LSGUF). Finally, we consider the case in which both the inner products and the scaling 
are chosen to minimize the least-squares error between the original frame and the resulting tight 
frame. The optimizing frame is the least-squares GU frame (LSGUF). 

In Section ^, we consider distance properties of GU frames, which may be of interest when using 
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GU frames for code design (group codes) |T^, In particular, we introduce a class of GU frames 
with strictly positive distance spectra for all choices of generating vectors. Such GU frames are 



shown to be generated by fixed-point-free groups [133 1. 

Before proceeding to the detailed development, in Section ^ we provide a brief introduction to 
frame expansions. 



2 Frames 

Frames, which are generalizations of bases, were introduced in the context of nonharmonic Fourier 
series by Duffin and Schaeffer Q (see also [^). Recently, the theory of frames has been expanded 
[^, 1^, in part due to the utility of frames in analyzing wavelet decompositions. 

Let {4>i,l < i < n} denote a set of n complex vectors in an m-dimensional Hilbert space TL. 
The vectors 0j form a frame for 7i if there exist constants A > and B < oo such that 

n 

A\\x\\^<^\{x,^i)\^<B\\x\\\ (1) 

i=l 

for all X G W In this paper, we restrict our attention to the case where m and n are finite. The 
lower bound in (||) ensures that the vectors (pi span 7i; thus we must have n > m. Since n < cxd 
the right hand inequality of (|^) is always satisfied with B = Y17=i 1 1*^*11^1 so that any finite set of 
vectors that spans W is a frame for 7i. In particular, any basis for 7^ is a frame for 7i. However, 
in contrast to basis vectors which are linearly independent, frame vectors with n > m are linearly 
dependent. If the bounds ^ = S in (|l|), then the frame is called a tight frame. If in addition 
A = B = 1, then the frame is called a normalized tight frame. The redundancy of the frame is 
defined as r = n/m, i.e., n vectors in an m-dimensional space. 
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The frame operator corresponding to the frame vectors {(pi, 1 < i < n} is defined as |17] 



s = Y,<l>^^^ = ^^*^ (2) 



1=1 



where $ is the matrix of columns (pi, and (•)* denotes the Hermitian transpose. Using the frame 
operator, (|^) can be rewritten as 

AWxW^ < {Sx,x) < BWxW^. (3) 

From it follows that the tightest possible frame bounds A and B are given hy A = miuj Aj(S') 
and B = maxjAj(S'), where {Xi{S),l < i < m} are the eigenvalues of the frame operator S. 
Throughout the paper, when referring to "the frame bounds" we implicitly assume the tightest 
possible frame bounds unless otherwise stated. 

If the vectors {(pi, 1 < i < n} form a frame for TC, then any x € Ti. can be expressed as a linear 
combination of these vectors: x = XlILi - n > m the coefficients in this expansion are not 
unique. A possible choice isQ Oj = {(pi,x) where (pi are the dual frame vectors [17] of the frame 
vectors (pi, and are given by 

^i = S~^(P^. (4) 
We can express directly in terms of 5 as Q 



A + B 



^T.{^-^S], (5) 



where A and B are the frame bounds of {</>«}. The choice of coefficients Oj = {(pi,x) has the property 
that among all possible coefficients it has the minimal /2-iiorm 22 1. 



There are other choices of dual frame vectors yi such that for any x £ TC, x = 



^We use the notation {x,y) = x*y. 



Specifically, with Y denoting the matrix of columns yi, any other choice corresponds to Y of the 



form H 



y = 5-i$ + T(/-$*5-^$) , (6) 

where T is an arbitrary matrix with bounded elements. However, the particular choice yi = (pi has 
some desirable properties. Besides resulting in the minimal /2-iiorm coefficients, in many cases the 
choice yi = (pi yields frame vectors that share the same symmetries as the original frame vectors. 
Specifically, in Section ^ we show that the dual frame vectors associated with a geometrically 
uniform (GU) frame are also GU, and in Section ^ we show that the dual frame vectors associated 
with a compound GU (CGU) frame are also CGU. Finally, in the case of a tight frame the dual 
frame vectors lead to a particularly simple expansion. Specifically, in this case S = Aim so that 
= {l/A)Im, and the dual frame vectors are {(j)i = {1/A)(pi,l < i < n}. Since a tight frame 



expansion of a signal is very simple, it is popular in many applications |17|. 

Suppose we are given a set of vectors {(pi, I < i < n} that form a frame for TC, with frame 
bounds A ^ B. It may then be desirable to construct a tight frame from these vectors. A popular 
tight frame construction is the canonical tight frame [ p!7| , |^, 26, |2^], first proposed in the 



context of wavelets in |^. The canonical tight frame vectors {^i, 1 < i < n} associated with the 
vectors {c/'j, 1 < i < n} are given by 

fi, = S''/^(f>i., (7) 
where 5"^/^ is the positive-definite square root of S~^. We can express S~^^'^ directly in terms of 



S as 



where A and B are the frame bounds associated with the frame {</>«}. Note, that {US~^^'^(j)i} with 
U an arbitrary unitary matrix yields a tight frame as well. The canonical tight frame, however, 
has the property that it is the closest normalized tight frame to the vectors {(pi} in a least-squares 



sense ^ . 

From (^) and (j^ we see that in order to compute the dual frame vectors and the canonical tight 
frame vectors associated with a frame {4>i}, we need to compute the matrices S^^ and S~^^^ and 
then apply them to each of the frame vectors (pi. In the next section, we introduce a class of frames 
that have strong symmetry properties called geometrically uniform (GU) frames. As we show in 
Section ^ the dual frame vectors and the canonical tight frame vectors associated with a GU frame 
are generated by a single generating function, and can therefore be computed very efficiently. 



3 Geometrically Uniform Frames 



A set of vectors S = {</>j, 1 < i < n} is geometrically uniform (GU) p9| , 18, 25] if every vector 
in the set has the form (pi = UiCp, where (p is an arbitrary generating vector and the matrices 
{Ui,l < i < n} are unitary and form an abelian group|^ Q. For concreteness we assume that 
Ui = I so that (pi = (p. The group Q will be called the generating group of S. 

Alternatively, a vector set is GU if given any two vectors (pi and (pj in the set, there is an 
isometry (a norm-preserving linear transformation) Zij that transforms (pi into (pj while leaving the 



set invariant Thus for every i, ZijCpi = (pj. Intuitively, a vector set is GU if it "looks the same" 



geometrically from any of the points in the set. Some examples of GU vector sets are considered 



m 



0. 



A set of vectors {(pi £ 7i,l < i < n} forms a geometrically uniform frame for TC, if the vectors 
{(pi} are GU and span TC. 

As we show in the following proposition, the frame bounds of a GU frame can be bounded by 
the norm of the generating vector. 

Proposition 1. Let S = {(pi = Ui(p, Ui E Q} be a geometrically uniform frame with frame bounds 



^That is, Q contains the identity matrix I; if Q contains Ui, then it also contains its inverse ^; the product 
UiUj of any two elements of Q is again in Q; and UiUj = UjUi for any two elements in Q |^ . 



A and B, where (j) is an arbitrary generating vector. Then A < < B. If in addition the 

frame is tight, then A = B = . 

Proof. We can express the frame operator corresponding to the frame vectors 0, as 

n 

S = J2Ui<l>^*U*. (9) 



i=l 



Then, 



Tr{S)=J2'^iUi^^*U:)=nM\^ (10) 



i=l 



so that 



J2HS)=Tr{S) = nMf. (11) 



1=1 



Therefore, 



-. m 

A = mm\i{S) < -y2HS) = -M\\ (12) 
1=1 



and 



^ m 

B = maxA,(5) > - V Ai(5) = (13) 
m ^—^ m 



□ 

Since = Uf^, the inner product of two vectors in 5 is 

= rU-^Uj<l> = s{U-^Uj), (14) 

where s is the function on Q defined by 

s{Ui) = <t>*Ui<t>. (15) 
For fixed i, the set U^^Q = {Uf^Uj,Uj £ Q} is just a permutation of Q since Uf^Uj G Q for 



all i,j |36| |. Therefore, the n numbers {s{Uf^Uj), 1 < j < n} are a permutation of the numbers 
{s{Ui), 1 < i < n}. The same is true for fixed j. Consequently, every row and column of the n x n 
Gram matrix G = {{(pi, (f)j)} is a permutation of the numbers {oj = s{Ui), 1 < i < n}. 

A matrix G whose rows (columns) are a permutation of the first row (column) will be called a 
permuted matriJ^ Thus, we have shown that the Gram matrix of a GU vector set is a permuted 
matrix. Furthermore, if the Gram matrix G = {{(f>i,4>j)} is a permuted matrix and in addition 
G = G^, then the vectors {(pi} are GU |37]. We therefore have the following proposition. 



Proposition 2. The Gram matrix G = {{(pi , (pj)} corresponding to a geometrically uniform vector 
set S = {(pi £ 7i,l < i < n} is a permuted matrix. Conversely, if the Gram matrix G = {{(pi, (pj)} 
is a permuted matrix, and {(pi,(pj) = {(pj,(pi) for all i,j, then the vectors {(pi} are geometrically 
uniform. If in addition the vectors {(pi} span TL, then they form a geometrically uniform frame for 

n. 

It will be convenient to replace the multiplicative group Q by an additive group Q to which 
Q is isomorphic^. Specifically, it is well known (see e.g., [Q) that every finite abelian group Q is 
isomorphic to a direct product Q of a finite number of cyclic groups: Q = Q = x • • • x Z„p, 
where Z„j is the cyclic additive group of integers modulo ut, and n = Yl^nt- Thus every element 
can be associated with an element q G Q of the form q — {qi, Q2, ■ ■ ■ , Qp)i where qt G ^nt! 
this correspondence is denoted by Ui <-> q. 

Each vector (pi = Ui(p is then denoted as (p{q), where Ui ^ q. The zero element = (0, 0, . . . , 0) G 
Q corresponds to the identity matrix I ^ Q, and an additive inverse —q S Q corresponds to a 



An example of a permuted matrix is 

ai a2 as 04 

a2 ai a4 as 

as 04 fli a2 

ai as a2 ai 



(16) 



* Two groups Q and Q' are isomorphic, denoted by Q = Q', if there is a bijection (one-to-one and onto map) 
V • Q ^ Q' which satisfies ip{xy) — (p{x)(p{y) for all x,y £ Q |3h]. 
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multiplicative inverse U- ^ = E Q. The Gram matrix is then the n x n matrix 



G = {W),<t>{q)),q', q&Q} = {s{q - q'),q', q £ Q}, (17) 
with row and column indices q',q E Q, where s is now the function on Q defined by 

s{q) = {m,<Piq))- (18) 

The Fourier transform (FT) of a complex-valued function ip : Q ^ C defined on Q = Z„j x 
• • • X Z„ is the complex- valued function ip : Q C defined by 



^(/t) = ^5^(/i,g)^(g), (19) 

where the Fourier kernel {h,q) is 

p 

= (20) 
t=i 

Here ht and qt are the kth components of h and q respectively, and the product htqt is taken as an 
ordinary integer modulo nj. 

The FT matrix over Q is defined as the n x n matrix = {^{h,q),h,q € Q}. The FT of a 
column vector ip = {<p{q),q € Q} is then the column vector (p = {(p{h),h € Q} given hy (p = Tip. 
Since T is unitary, we obtain the inverse FT formula 



p = r'p = < 



-^Y.{h,qYp{h),q^Q)- 

heQ I 



(21) 
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As we show in the following theorem, the FT matrix plays an important role in defining GU 
frames. 

Theorem 1. A set of vectors {(pi, 1 < i < n} in an m- dimensional Hilbert space 7i is geometrically 
uniform if and only if the Gram matrix G = {{(pi, (pj)} is diagonalized by a Fourier transform matrix 
T over a finite product of cyclic groups Q. The vectors form a geometrically uniform frame 
for TC if in addition G has rank m. 

Proof. The vectors {(pi} form a frame for 7i if and only if they span 7i, which implies that the rank 
of G must be equal to m. 

For a GU vector set with generating group Q = Q, the FT over Q diagonalizes the Gram matrix 
G ||25|. Thus, to complete the proof of the theorem we need to prove that if G is diagonalized by 
a FT matrix over the group Q, then the vector set {(j){q),q G Q} is GU. 

Let ^ be the matrix of columns (j){q), so that G = Since diagonalizes G, G has an 

eigendecomposition of the form G = TDT* for a diagonal matrix D with diagonal elements di, 
where the first m diagonal elements may be non-zero and the remaining diagonal elements are all 
zero. Then $ has an SVD of the form $ = UTiT* where U is an arbitrary unitary matrix and S is 
an m X n diagonal matrix with diagonal elements {o", = ^/dl, 1 < i < m}. 

Let {f{q),q G Q} denote the columns of J^* . From the definition of the components of /(O) 
are all equal 1/y/n, and f{q) = B{q)f{0) where B{q) is a diagonal unitary matrix with diagonal 
elements {{h,q),h G Q}, where {h,q) is given by (pc|). Then, 

^{q) = U^f{q) = UJ:Biq)fiO) = C/B(g)E/(0) = UB{q)U*(t>, (22) 

where (p = C/S/(0), and where we used the fact that diagonal matrices commute. If we now define 
U{q) = U B{q)U* , then we have that (p[q) = U{q)(p where the matrices {U{q),q G Q} are unitary. 
We now show that the group Q = {U{q),q G Q} is an abelian group. First, we have that 
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C/(0) = UB{0)U* = UU* = I so that I e Q. Next, U-^{q) = UB-^{q)U* = UB{-q)U* so that 
U'^{q) e Q since -q E Q. Finally, U{q)U{h) = UB{q)B{h)U* = UB{h)B{q)U* = U{h)U{q) 
since diagonal matrices commute, and U{q)U{h) G Q since B{h)B{q) = B(h + q) = B{q') for some 

q'eQ. 

We therefore conclude that (j)(q) = U{q)(j) where the matrices U{q) are unitary and form an 
abelian group, so that the vectors (j){q) are geometrically uniform. □ 

As a consequence of Theorem |^ we have the following corollary. 

Corollary 1. A set of vectors {(pi £ 7i,l < i < n} is geometrically uniform if and only if the 
matrix ^ of columns (pi has an SVD of the form $ = UTiT* , where U is an arbitrary unitary 
matrix, E is an arbitrary diagonal matrix with diagonal elements {ai, 1 < i < m}, and T is a FT 
matrix over a direct product of cyclic groups. In addition, the vectors {0j} form a geometrically 
uniform frame for 7i if they span 7i or, equivalently, if ai ^ for 1 < i < m. 



4 Dual and Canonical Tight Frames Associated With GU Frames 



In Section 4.1, we show that the dual frame vectors and the canonical tight frame vectors associated 
with a GU frame are also GU. This property can then be used to compute the dual and canonical 
tight frames very efficiently. Further properties of the canonical tight frame vectors are discussed 



in Section 4.2 



4.1 Constructing the dual and canonical tight frames 

Let S = {(pi = Ui(p, Ui € Q} be a GU frame generated by a finite abelian group Q of unitary 
matrices, where (p is an arbitrary generating vector. Then the frame operator S defined by (^) 
commutes with each of the unitary matrices Ui in the generating group Q. Indeed, expressing the 
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frame operator as 



S = Y,U^c|>c^*U*, (23) 



i=l 



we have that for all j, 

n 

i=l 

n 

= u,^u*u,(t>ru:uj 

i=l 
n 

i=l 

= UjS, (24) 

since {U*Ui, 1 < i < n} is just a permutation of Q. 

If S commutes with Uj, then from (|5|) and (^) we have that and S"^^/^ also commute with 
Uj for all j. Thus, 

= s-^4>^ = s~^Ui(j) = ms-^4> = Ui4>, (25) 

where (j) = S~^<j), which shows that the dual frame vectors {(^j = S~^(f)i] are GU with generating 
group equal to Q. 
Similarly, 

= 5-1/2 ^ = C/,S-i/20 = u.^^ (26) 

where // = S~^^'^4>, which shows that the canonical tight frame vectors {fii = S~^^'^(pi} are also GU 
with generating group Q. 

Therefore, to compute the dual frame vectors or the canonical tight frame vectors all we need 
is to compute the generating vectors (j) ^-nd /x, respectively. The remaining frame vectors are then 
obtained by applying the group Q to the corresponding generating vectors. 

We now show that the generating vectors can be computed very efficiently using the FT. From 
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Corollary |^ we have that $ has an SVD of the form 



$ = f/S^* = (T{h)u{h)T* (h) . (27) 

Here S is a diagonal matrix with diagonal elements = n^^^y^s{h), h G Q} where {s{h), h G Q} 

is the FT of {s{q),q £ Q}, U is the matrix of columns u{h), where 



u{h) 



^T{h)la(h) = Mh)/a{h), if aih) ^ 0; 

(28) 

0, otherwise 



with 

m = ^y^{Kq)<t^{q) (29) 

denoting the hth element of the FT of $ regarded as a row vector of column vectors, $ = {(l){q),q G G} , 
andT* ={^{h,q)*,h,qeQ} has rows J^*{h) = {^{h,q)* ,q e Q}. 
It then follows that 

^ h&x ^ ' 

where /i G T if cj(/i) 7^ 0. Similarly, 



5-1/2^ = ^ J] n(/i). (31) 



In 

hex 
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We summarize our results in the following theorem: 

Theorem 2 (GU frames). Let S = {(pi = Uicjy^Ui £ Q} be a GU frame generated by a finite 
abelian group Q of unitary matrices, where (p is an arbitrary generating vector, and let 4? be the 
matrix of columns (pi. Let Q be an additive abelian group isomorphic to Q, let {(p{q),q £ Q} be the 
elements of S under this isomorphism, and let T be the Fourier transform matrix over G. Then 

1. the dual frame vectors {i^j, 1 < i < n} are geometrically uniform with generating group Q and 
generating vector 4> = (1/a/^) where 

(a) {cr{h) = n^^^^/s{h), h G Q} are the singular values of ^, 

(b) {s{h), h £ Q} is the Fourier transform of the inner-product sequence {(0(0), (piq)), q G Q}, 

(c) Z is the set of indices h £ Q for which a{h) / 0, 

(d) u{h) = ^{h)/a{h) for h£l, 

(e) {(j){h),h £ Q} is the Fourier transform of{(j){q),q £ Q}, 

2. the canonical tight frame vectors {/Xj, 1 < i < n} are GU with generating group Q and 
generating vector fi = {l/^/n) J2hel'^(^)' 

3. the frame bounds of the frame {^i,! < i < ^^} are given by A = y/nmiuhi^x s{h) and B = 
/nmax/igj s{h). 



An important special case of Theorem g is the case in which the generating group Q is cyclic 
so that Ui = V^^^, 1 < i < n, where y is a unitary matrix with = /. A cyclic group generates 
a cyclic vector set S = {(pi = V^~^(p, 1 < i < n}, where (p is arbitrary. If Q is cyclic, then G is 
a circulant matrix]^, and Q is the cyclic group The FT kernel is then {h, g) = g-^'^^^f/" for 



circulant matrix is a matrix where every row (or column) is obtained by a right circular shift (by one position) 

ao a2 ai 
a\ ao a2 
a-z ai ao 



of the previous row (or column). An example is: 
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h,g £ Z„, and the FT matrix reduces to the n x n DFT matrix. The singular values of $ are 
then n^/"^ times the square roots of the DFT values of the inner products {{(pi, (f>j)A ^ J ^ n}. 

4.2 Properties of the canonical tight frame 

The canonical tight frame vectors fii corresponding to the frame vectors (pi have the property 
that they are the closest normalized tight frame vectors to the vectors in a least-squares sense 
11^. Thus, the vectors /Xj are the normalized tight frame vectors that minimize the least- 
squares error 

n 

^{(pi- fii,(l)i- Hi). (32) 
1=1 

We now show that when the original frame vectors (pi are GU with generating group Q, the 
canonical tight frame vectors have the additional property that from all normalized tight frame 
vectors they maximize 

n 

R^^. = Y.m,^^i)\^■ (33) 

i=l 

Maximizing i?,^^ may be of interest in various applications. For example, in a matched-filter 
detection problem considered in |jl2|, R^^ represents the total output signal-to-noise ratio. As 
another example, in a multiuser detection problem considered in |13|, maximizing R^^ has the 



effect of minimizing the multiple-access interference at the input to the proposed detector. 

To obtain a more convenient expression for R^^, let <I> and M denote the matrices of columns 
(pi and fii, respectively. Since the vectors /ij form a normalized tight frame for 7i, M satisfies 



MM* = Im. (34) 

From Corollary 0, $ has an SVD of the form <I> = UT,J-* , where U is unitary, is the FT matrix 
over the additive group Q to which Q is isomorphic, and S is an m x n diagonal matrix with 
diagonal elements cTj > 0. From (]3^ ) it follows that M can be written as M = UIZ* where Z is an 
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arbitrary unitary matrix and / is an m x n diagonal matrix with diagonal elements all equal to 1. 
Let fi and Zj denote the columns of J-'* and Z* , respectively. Then we can express R^f^^ as 

n n n 

Rc,, = E = E l^^**^- f^>^)l' = ^\{f^,^^^)\^ (35) 

i=l i=l 1=1 

where S is an n x n diagonal matrix with the first m diagonal elements equal to CTj, and the 
remaining diagonal elements are all equal to 0. 

Our problem then reduces to finding a set of orthonormal vectors Zi that maximize | {fi, Szj) p, 
where the vectors fi are also orthonormal. Using the Cauchy-Schwarz inequality we have that 

n n 

R^^ = J] < S/,)(z„ Sz,), (36) 

i=l i=l 

with equality if and only if T,^^'^ fi = CiT}^"^ Zi for some Cj. In particular, we have equality for Zi = fi. 
Since the components of the vectors fi all have equal magnitude 1/-Jn, {fi, T,fi) = X^fcLi (^k=oi 
for all i, and ( p6| ) reduces to 

n 

R<Pf, < a^(zi,Szi) = aTr(Z*SZ) = aTr(S) = na^ (37) 

i=l 

with equality if fi = Zi. 

The normalized tight frame vectors that maximize i?,^^ are then the columns of M = U IJ-* = 
5-i/2(j) where S = <!>$*, and are equal to the canonical tight frame vectors. 
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5 Example Of A GU Frame 

We now consider an example demonstrating the ideas of the previous section. 

Consider the frame vectors 0i = l/2[^/3 - l]*,(j)2 = 1/2[V3 l]*,(l)3 = l/2[-\/3 1]*,04 = 
l/2[-y/3 - 1]*, depicted in Fig. |l|. 

The corresponding Gram matrix is given by 



1 


0.5 


-1 


-0.5 


0.5 


1 


-0.5 


-1 


-1 


-0.5 


1 


0.5 


-0.5 


-1 


0.5 


1 



(38) 



which is a permuted matrix with G = . From Proposition ^ it follows that the vectors (pi are 
GU. Since the vectors (pi also span TZ'^, these vectors form a GU frame for TZ^. 




Figure 1: Example of a GU frame. 



The vectors (pi can be expressed as {(pi = Ui(p), 1 < i < 4}, where (p = (pi and the matrices 
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{Ui, 1 < i < 4} are unitary, form an abelian group Q, and are given by 





1 




-1 




-1 


Ul=l2, U2 = 














-1 




-1 




1 



The multiplication table of the group Q is 



If we define the correspondence 





Ul 


U2 


Us 


Ua 


Ui 


Ui 


U2 


Us 


U4 


U2 


U2 


Ul 


Ua 


Us 


Us 


Us 


Ui 


Ul 


U2 


f/4 


Ui 


Us 


U2 


Ul. 



(39) 



(40) 



C/i^(0,0), t/2^(0,l), Us ^{1,0), [/4^(1,1), 
then this table becomes the addition table of Q = 7^2 x Z2: 





(0,0) 


(0,1) 


(1,0) 


(1,1) 


(0,0) 


(0,0) 


(0,1) 


(1,0) 


(1,1) 


(0,1) 


(0,1) 


(0,0) 


(1,1) 


(1,0) 


(1,0) 


(1,0) 


(1,1) 


(0,0) 


(0,1) 


(1,1) 


(1,1) 


(1,0) 


(0,1) 


(0,0). 



(41) 



(42) 



Only the way in which the elements are labeled distinguishes the table of (El) from the table of 
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(^0[); thus Q is isomorphic to Q. Over Q = TLi x Z2, the FT matrix T is the Hadamard matrix 



1 1 

-1 1 

1 -1 

1 -1 



(43) 



From Theorem |2| the dual frame vectors and the canonical tight frame vectors are also GU with 
generators and respectively, whose equations are given in the theorem. Thus, to compute the 
dual and canonical tight frame vectors we compute these generators and then apply the group Q. 

We first determine the FT of the first row of G denoted by s: 



s = J^s = -[0 3 1]*. 



(44) 



Using Theorem g, it follows from ( |4^ ) that the frame bounds are given by ^ = 1 and 5 = 3. Next, 
we compute the vectors which are the columns of 



1.7 

-1 



(45) 



Using the equations of the theorem, we then have that = [0.3 — 0.5]* and = [0.5 — 0.5]*. By 
applying the group Q to these generators we obtain that the dual frame vectors are the columns of 



0.3 0.3 -0.3 -0.3 
-0.5 0.5 0.5 -0.5 



(46) 
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and the canonical tight frame vectors are the columns of 



M 



0.5 0.5 -0.5 -0.5 
-0.5 0.5 0.5 -0.5 



(47) 



Comparing ( [dq) and ( |47| ) with the original frame vectors 4>i, it is evident that the dual and 
canonical tight frame vectors have the same symmetries as the original frame vectors, as illustrated 
in Fig. |. 




Figure 2: Symmetry property of the frame vectors (pi, the dual frame vectors vectors (pi, and the 
canonical tight frame vectors /Xj. (pi are the columns of $ given by (p6|), and /ij are the columns of 
M given by (p7|). The frame vectors, dual frame vectors, and the canonical tight frame vectors all 
have the same symmetry properties. 
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6 Compound GU Frames 

In Section ^ we showed that the dual and canonical tight frame vectors associated with a GU frame 
are themselves GU and can therefore be computed using a single generator. In this section, we 
consider a class of frames which consist of subsets that are GU, and are therefore referred to as 
compound geometrically uniform (CGU) frames. As we show, the dual and canonical tight frame 
vectors associated with a CGU frame share the same symmetries as the original frame and can be 
computed using a set of generators. 

A set of frame vectors {(j)ik, 1 < i < /, 1 < < r} is CGU if (pik = Uicpk for some generating 
vectors {(pk^ ^ ^ k < r}, and the matrices {Ui, 1 < i < 1} are unitary and form an abelian group Q. 

A CGU frame is in general not GU. However, for every k, the vectors {(f)ik, 1 < « < ^ GU 
vector set with generating group Q. 

A special case of CGU frames are filter bank frames studied in [^, ^, in which Q is the group 
of translations by integer multiples of the subsampling factor, and the generating vectors are the 
filter bank synthesis filters. We may therefore view CGU frames as an extension of filter bank 
frames to the more general group case. 

As we show in the following proposition, the frame bounds of a CGU frame can be bounded by 
the sum of the norms of the generating vectors. 

Proposition 3. Let S = {(pn. = Ui(pk, be a compound geometrically uniform 

frame with frame bounds A and B, where {(pk, 1 < k < r} is an arbitrary set of generating vectors. 
Then A < :f^Yl^k=i \ \4'k\\'^ ^ B. If in addition the frame is tight, then A = ;^ X]fc=i ll'/'fclP- 

Proof. We can express the frame operator corresponding to the frame vectors (pi as 

I r 

s = Y,Y.^,<Pkmi (48) 

i=l k=l 
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Then, 



I r 



Tr(5) = ^^TV([/,</-,<^^i7*) = i^||<^fc||2, (49) 



i=l k=l k=l 



SO that 



J2HS) = Tr{S) = lY,m\^ (50) 



i=l k=l 



Therefore, 



and 



m J r 

A = mmX,{S) < - VA,(5) = - (51) 



m — ' m 

1=1 fe=i 



1 "' I ^ 

B = max Ai(5) > - V Xi{S) = - V (52) 



m — ' m 

1=1 k=l 



□ 



6.1 Dual and canonical tight frames associated with CGU frames 

We now show that the dual and canonical tight frames associated with a CGU frame are also CGU. 
Expressing the frame operator as 

I r I / r \ 



i=l k=l i=l \k=l / 

for all j we have that, 

i=l \k=l J 

i=l \k=l J 

i=l \k=l J 
= UjS, (54) 
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since {UjUi, 1 < i < 1} is just a permutation of Q. Thus S commutes with Uj, so that S ^ and 
also commute with Uj for all j. Then, the dual frame vectors (pi^ of the vectors (pi^ are given 

by 

4>ik = S~^(f)ik = S~^Ui(j)k = UiS^'^^k = Ui4>k, (55) 

where 4>k = S~^4>k, which shows that the dual frame vectors = S~^<pih} are CGU with 

generating group equal to Q. 
Similarly, 

fiik = S-^/^ct>ik = S-^/^Uicl>k = UiS-^/^ct>k = U,^k, (56) 

where Hk = S~^f'^4>k, which shows that the canonical tight frame vectors {yiik = S~^/'^(pik} are also 
CGU with generating group Q. 

Therefore, to compute the dual frame vectors or the canonical tight frame vectors all we need 
is to compute the generating vectors {<pkA < k < r} and {^fc, 1 < k < r}, respectively. The 
remaining frame vectors are then obtained by applying the group Q to the corresponding set of 
generating vectors. 

6.2 CGU frames with GU generators 

A special class of CGU frames is CGU frames with GU generators in which the generating vectors 
{^fc) 1 < A; < r} are themselves GU. Specifically, {(pk = V^cp} for some generator cp, where the 
matrices {V^, 1 < A; < r} are unitary, and form an abelian group Q. 

Now suppose that Up and Vt commute up to a phase factor for all t and p so that UpVt = 
VtUpC^^^P'^^ where 9 {p,t) is an arbitrary phase function that may depend on the indices p and t. In 
this case we say that Q and Q commute up to a phase factor. Then for all p, t, 

SUpVt = ^Ui(j2^kCpcp*vAu*UpVt 

1=1 Vfe=i / 
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i=l \k=l / 

i=l \k=l / 

UpVtY^^ (Yyt*yk<t^<t^*ykyi] u^ 



i=i \k=i 

I / r 



UpVtYu^[Y^k't^'t'*^k]u: 
i=i \k=i J 

UpVtS. (57) 



The dual frame vectors (pn. of the vectors (pn. are then given by 

^ik = S-^<Pik = S~^U,Vk(P = UiVkS-'cP = UiVkl (58) 

where (j) = S~^(j). Similarly, 

l^ik = S-^'^ik = S-^'^UiVk4> = UiVkS-^/^cf) = UiVkii (59) 

where ^ = S~^^'^4>. Thus even though the frame is not in general GU, the dual and canonical tight 
frame vectors can be computed using a single generating vector. 

As we now show, in the special case in which ^ = so that UiVk = VkUi for all i, k, the resulting 
frame is GU. To this end we need to show that the unitary matrices Q' = {Qik = UiV^} form an 
abelian group. First, QikQjt = UiVkUjVt = UiUjVkVt- Since UiUj G Q and VkVt G Q, QikQjt G Q' ■ 
Next, Q;,} = V^U* = U*V^ G Q' since U* e Q and G Q. Also, / G Q' since I G Q and / G C?. 
Finally, QikQjt = UiVkUjVt = UiUjVtVk = UjVtUiVk = QjtQik- 

A special case of CGU frames with GU generators for which Q and Q commute up to a phase 
factor are Weyl-Heisenberg (WH) frames [17, ^, If the WH frame is critically sampled, then 



9{p, t) = and the WH frame reduces to a GU frame. In the more general oversampled case, 
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^iP^t) 7^ 0- CGU frames with GU generators can hence be viewed as a generahzation of WH 
frames to the group case. 

To summarize, we have the following theorem: 

Theorem 3 (Compound GU frames). Let S = {(pik = Ui<pkA < i < 1,1 < k < r} be a 

compound geometrically uniform frame with generating vectors {(f)k, 1 < A; < r} and generating 
group Q, and let S be the frame operator corresponding to the frame vectors {(pik}- Then 

1. the dual frame vectors {4>ik, I < i < 1,1 < k < r} are compound geometrically uniform with 
generating group Q and generating vectors {(pk = S~^4>k, 1 < k < r}; 

2. the canonical tight frame vectors {nik,^ < i < l,^ < k < r} are compound geometrically 
uniform with generating group Q and generating vectors {fi^ = 5"^/^^^, 1 < A; < r}. 

// in addition the generating vectors {0^ = Vkcp, 1 < k < r} are geometrically uniform with UiV^ = 
VkUiC^^^''''^ for all i, k, then 

1- 4>ik = UiVk4> where 4> = S~^(f); 

2- fJ-ik = UiVkH where fi = <S'~^/^0; 

3. if in addition 6{i,k) = for all i,k, then the vectors < i < l,^ < k < r} form a 
geometrically uniform frame. 
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7 Pruning GU Frames 

In applications it is often desirable to know how a frame behaves when one or more frame elements 
are removed. In particular, it is important to know or to be able to estimate the frame bounds of 
the reduced frame. In general, if no structural constraints are imposed on a frame this behavior will 
depend critically on the particular frame elements removed. For example, removal of a particular 
frame element may destroy the frame property so that the remaining vectors do not constitute a 
frame anymore, whereas if a different element is removed the remaining vectors may still constitute 
a frame. 

One of the prime applications of frames is signal analysis and synthesis, where a signal is ex- 
panded by computing the inner products of the signal with the frame elements. The resulting 
coefficients are subsequently stored, transmitted, quantized or manipulated in some way. In par- 
ticular, a coefficient may be lost {e.g., due to a transmission error) which results in a reconstructed 
signal that is equivalent to an expansion using a pruned frame obtained by removing the corre- 
sponding frame vector. 

Recently, there has been increased interest in using frames for multiple-description source coding 
where a signal is expanded into a redundant set of functions and the resulting coefficients are 
transmitted over a lossy packet network where one or more of the coefficients can be lost because 



a packet is dropped [14|. The goal of multiple description source coding is to ensure a gradually 
behaving reconstruction quality as a function of the number of dropped packets. When using frames 
in this context, the reconstruction quality is often governed by the frame bound ratio of the pruned 
frame. If the packets are dropped with equal probability, then it is desirable that the frame bound 
ratio should deteriorate uniformly irrespectively of the particular frame element that is removed. In 
the following, we show that GU frames have this property. We will furthermore demonstrate that 
if the original frame is a tight GU frame, then the frame bound ratio of the pruned frame obtained 
by removing one frame element can be computed exactly. We also consider the case where sets of 
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frame elements are removed. 

Theorem 4 (Pruned GU frames). Let S = {^^ = Ui(t),Ui G Q} 6e a geometrically uniform 
frame generated by a finite abelian group Q of unitary matrices, where (j) is an arbitrary generating 
vector. Let $ be the matrix of columns (pi, and let S = be the corresponding frame operator. 
Let S{j) = {(pi = Ui(p,Ui E Q,i ^ j} be the pruned set obtained by removing the element (pj. 

Then the eigenvalues of the frame operator corresponding to the pruned set do not depend on the 
particular element (pj removed. 

Proof. The frame operator corresponding to the pruned frame is given by 

n 

S{j) = E Uicl>cP*U* - Uj(P(P*U;. (60) 
1=1 

Since Uj is unitary, the eigenvalues of S{j) are equal to the eigenvalues of UjS{j)Uj. But, 

n n 

U*S{j)Uj = U* J2 Ui^^*U*Uj -(P(P* = J2 Ui(P(P*U* -(P(P* = S- (61) 

i=l 1=1 
Since U*S{j)Uj is independent of j, the eigenvalues of S{j) do not depend on j. □ 

In general it is difficult to provide estimates on the frame bounds of the pruned frame. However, 
in the special case where the original GU frame is a tight frame, these bounds can be determined 
exactly. 

Corollary 2 (Pruned Tight GU frames). Let S = {(pi = Ui(p,Ui G Q} 6e a geometrically 
uniform tight frame generated by a finite abelian group Q of unitary matrices, where (p is a unit 
norm generating vector. Let $ be the matrix of columns (pi, and let S = be the corresponding 
frame operator. Let S{j) = {(pi = Ui(p,Ui E Q,i ^ j} be the pruned set obtained by removing the 
element (pj . Then the eigenvalues of the frame operator corresponding to the pruned set are given 
by Xi = ^ — 1 and Aj = ^ , 2 < i < n, independent of (pj . 
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Proof. Since 5 is a tight frame with ||i^|| = 1, from Proposition |l] the frame bound A = n/m and 
S = {n/m)Im- Then, 

7? 

S{j) = -Im - Uj<l)ct>*U*, (62) 

and 

U;S{j)U, = -!„,- ^^*. (63) 

Since = 1, has one eigenvalue equal to 1, and the remaining eigenvalues equal to 0. The 
eigenvalues of S{j) are therefore given by Ai = ^ — 1 and Xi = 2 < i < n. □ 

An immediate consequence of Corollary Q is that the frame bound ratio of the pruned frame is 
given by B/A = 1/(1 — m/n), which is close to 1 for large redundancy r = n/m. 
We next consider the case where multiple frame elements are removed. 

Corollary 3. Let S = {(pi = Uicj), Ui G Q} be a geometrically uniform frame generated by a finite 
abelian group Q of unitary matrices, let $ be the matrix of columns (pi, and let S = <I><5* be the 
corresponding frame operator. Let J' be a set of indices, and let J{k) denote the set of indices 
i such that Ui = U^Uj for fixed k and j £ J . Let S{k) = {(pi = Ui(p,Ui £ Q,i ^ J^ik)} be a 
pruned set obtained by removing the elements (pi with i G J{k). Then the eigenvalues of the frame 
operator corresponding to the pruned set are independent of k. 

Proof. The frame operator corresponding to the pruned frame is given by 



s{k) = yuM*u* -Uk 



1=1 



UI (64) 



Then, 



UlS{k)Uk = Ul^Ui(P(P*U:Uk-Y,Uj(P(P*U* 

i=l j£j 
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n 

i=l j&J 

is independent of k, and consequently the eigenvalues of S{k) do not depend on k. □ 

To conclude this section, GU frames have strong symmetry properties in the sense that removing 
any one of the elements leads to a vector set with bounds independent of the particular element 
removed. Moreover, if the original frame is tight, then we can compute the bounds of the pruned 
frame exactly. 

8 Constructing GU Frames 

Suppose we are given a set of vectors {(pi, 1 < i < n}, that form a frame for an m-dimensional space 
TC. We would like to construct a geometrically uniform frame {0j,l < i < n} from the vectors 

Wi}- 

From Theorem ^ it follows that the vectors form a GU frame if and only if the Gram matrix 
G has rank m, and is diagonalized by a FT matrix over a finite product of cyclic groups. There 
are many ways to construct a frame from a given set of frame vectors <pi that satisfy these properties. 
For example, let F be the matrix of columns (pi, and let F have an SVD F = QAV* , where A is 
a diagonal matrix with diagonal elements Aj. Then the columns of <I> = FVEJ^* = QATiJ^* form 
a GU frame, where J- is any FT matrix over a product of cyclic groups, and S is an arbitrary 
diagonal matrix with diagonal elements CTj > 0. The frame bounds of the resulting GU frame are 
given hy A = miuj Xfaf and B = maxj Xfcrf, so that we can choose the diagonal matrix S to control 
these bounds. In particular, choosing S = / we have that the columns of <I> = FVJ^* = QKT* 
form a GU frame. This choice has the property that the frame bounds of the GU frame are equal 
to the frame bounds of the original frame. 

We now consider the problem of constructing an optimal GU frame. Specifically, let < 
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i < n} denote a frame for TC, and suppose we wish to construct a GU frame {(pi} from the vectors 
{(fi}- A reasonable approach is to find a set of vectors <j)i that span TC, and are "closest" to the 
vectors ipi in the least-squares sense. Thus we seek vectors (pi that minimize the least-squares error 
E, defined by 

n 

E = Y,{ei,e^), (65) 



i=l 



where e,- denotes the ith error vector 



ei = ipi-(pi, (66) 



subject to the constraint that the vectors (pi form a GU frame. 

If the vectors (pi are GU, then their Gram matrix G = is a permuted matrix with rank m, 
diagonalized by a FT matrix J^. Thus, the inner products {{(pi,(pj)} must satisfy 

{{(Pi, (Pj), l<j<n}= (3^Pi{a,, l<j< n}, (67) 

where Pi{aj, 1 < j < n} is a permutation of the numbers {aj, 1 < j < n}, /5 > is a scaling factor, 
and the numbers {aj,l < J < Jt-} are chosen such that the matrix R whose ith row is equal to 
1 < J < J^} is Hermitian, non-negative definite and diagonalized by a FT matrix J^. 

In our development of the optimal GU frame vectors we assume that the permutations Pi in ( |67| ) 
are specified. Since these permutations determine the additive group Q over which the FT matrix 

is defined, we assume that is specified. We then consider three different constraints on the 
vectors (pi. First we consider the case in which both the numbers {aj, 1 < j < n} and the scaling 
P in (^) are known. The GU frame minimizing the least-squares error E of (|65|)-(]66|) subject to 



this constraint is derived in Section B.l, and is referred to as the scaled-constrained least-squares 
GU frame (SC-LSGUF). Next, we consider the case in which the numbers {aj,l < j < n} in 
( |67|) are known, and the scaling /3 is chosen to minimize E. The resulting GU frame is referred 



to as the constrained least-squares GU frame (C-LSGUF), and is derived in Section B.2. Finally, 
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in Section |8.3| we consider the more general case in which both the numbers {aj, 1 < j < n} and 
the scaling (3 in ( |67| ) are chosen to minimize E. The resulting GU frame is referred to as the 
least-squares GU frame (LSGUF). 

8.1 Scaled-constrained least-squares GU frame 

We first consider the case in which the numbers {aj, 1 < j < n} and the scaling (3 in (|67|) are 
known. Thus, we seek the set of vectors {<pi} that minimize the least-squares error E of (|65|)-(|66|) 
subject to the constraint 

= = PiTAT*, (68) 

where $ is the matrix of columns /3o is a known scaling factor, and R is the matrix whose ith 
row is equal to Pi{aj, 1 < j < n} where the numbers {aj, 1 < j < n} are given such that R is 
diagonalized by JF, and A is a diagonal matrix with diagonal elements {aj = n^^'^aj, 1 < j < n} 
where {aj, 1 < i < n} is the FT of the sequence {uj, 1 < j < n}. From (|68[), the frame bounds of 
the vectors are given hy A = I3q miuj Uj and B = (3q maxj Uj. 

This problem has been solved in the context of general least-squares inner product shaping 
The optimal SC-LSGUF vectors are the columns of given by 

$ = /3oC/F*SJ^*, (69) 

where U and V are the right-hand unitary matrix and left-hand unitary matrix respectively in the 
SVD of FJ^T,*, F is the matrix of columns ipi, and S is an m x n diagonal matrix with diagonal 
elements y/oi for values of i for which ^ 0. 

If FRF* = PIFTAT*F is invertible, then we may express $ as 

S = j3Q{FJ'AJ'*F*)~^''^FJ'AF* = l3o{FRF*)-^/'^FR. (70) 
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8.2 Constrained least-squares GU frame 

We now consider the case in which the numbers {aj, 1 < j < n} are known, but the scahng P is 
not specified. Thus, we seek a set of vectors {(pi} that minimize the least-squares error E subject 
to 

= /?2^, (71) 

where (3 > 0. 

The least-squares error E of (|65|)-(|6^) may be expressed as 



E = Tr(($-F)*($-F)) 

= Tr(^>*$) + Tr(F*F) - 25R{Tr(F*^>)} 

= P^Tt{R) + Tt{F*F) -2^{Tr{F*^)}. (72) 
Let $ = (l//3)$. Then minimizing E is equivalent to minimizing 

E' = (3^Tr{R) - 2/?K{Tr(F*$)}, (73) 

subject to 

= R. (74) 
To determine the optimal matrix <I> we have to minimize E' with respect to /? and <I>. Fixing /? 



and minimizing with respect to ^, the optimal value of $ is given by the SC-LSGUF of Section 3.1 
with scaling = 1, so that 

$ = UV*^T*. (75) 
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If FRF* is invertible, then 

$ = {FJ'AT*F*)-^^^FrAT* = {FRF*)-^/'^FR. (76) 
Substituting <I> back into ([73|), and minimizing with respect to /?, the optimal value of j3 is given 



^ 5R{Tr(F*$)} Tr(C/F*SJ^*) 



by 

Ti{R) Tt{R) 



which in the case that FRF* is invertible reduces to 

Tr{{FRF*)y^) 

^- MR) • 

The C-LSGUF vectors are then the columns of <I> given by 

S = pUV*^J^*, (79) 
where /3 is given by ([77|). If FRF* is invertible, then 

$ = p{FRF*y^/'^FR, (80) 

where /3 is given by (|78|). 
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We summarize our results regarding constrained optimal GU frames in the following theorem: 

Theorem 5 (Constrained least-squares GU frames). Let {fi} be a set of n vectors in an m- 
dimensional Hilbert space Ti that span TL, and let F he the matrix of columns fi. Let {cpi} denote 
the optimal n GU frame vectors that minimize the least-squares error defined by subject 
to the constraint j^dlp, and let $ be the matrix of columns (f>i. Let R be the matrix with ith row 



equal to Pi{aj, 1 < j < n}, so that is the Gram matrix of inner products {{(pi, and let T 
be the Fourier transform matrix that diagonalizes R. Let A be the diagonal matrix with diagonal 
elements {aj = n^^'^aj, 1 < j < n} where {aj, 1 < j < n} is the FT of the sequence {aj, 1 < j < n}, 
let T, be an m X n diagonal matrix with diagonal elements yjai for values of i for which 7^ 0, 
and let U and V be the right-hand unitary matrix and left-hand unitary matrix respectively in the 
SVD ofFJ^T.*. Then, 

1. if f3 = Po is given, then ^ = PqUV*'SJ-'* . If in addition FRF* is invertible, then $ = 
(3o{FJ^AF*F*)-^/'^FFAJ^* = l3o{FRF*)-^''^FR; 

2. if (3 > is chosen to minimize E, then $ = /3UV*TiF* , where (3 is given by jT^j. // 
in addition FRF* is invertible, then $ = /3{FFAF* F*)-^/^FFAJ'* = P{FRF*)~^/^FR, 



where (3 is given by ^7q). 
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8.3 Least-squares GU frame 

We now consider the least-squares problem in which both the scaling factor (3 and the numbers 
{oj, 1 < j < f^-} in ( |67|) are chosen to minimize E. Thus, we seek a set of vectors that minimize 
the least-squares error E of (|65|)-(|66|) subject to 

{{ct)i, (t>j), 1 < j < n} = P,{aj, l<i< n}, (81) 

where Pj{aj, 1 < j < n} is a known permutation of the unknown numbers {oj, 1 < J < n}, chosen 
such that the matrix R whose iih. row is equal to Pi{aj, 1 < i < ?^} is Hermitian, non-negative 
definite and diagonalized by a FT matrix J^. 

This problem has also been considered in the context of general least-squares inner product 
shaping [^]. It can be shown that the solution involves solving a problem of the form 

m 

max ^ (82) 

i=l 

subject to 

{yi,yj) = Sij, (83) 

where the vectors {xi} are known and are a function of the given vectors {fi}- 

As we now show, this problem is equivalent to a quantum detection problem, for which there 
is no known analytical solution in the general case. 

8.3.1 Connection with quantum detection 

In a quantum detection problem, a system is prepared in one of m known (pure) states that are 
described by vectors Si in a Hilbert space 7i, and the problem is to detect the state prepared 
by performing a measurement on the system. The measurement is described in terms of a set of 



37 



orthogonal measurement vectors Qi. Given a set of measurement vectors qt, and assuming equal 
prior probabilities on the different states, the probability of detection is given by 



^ m 



(84) 



i=l 



Comparing (^) with (82) we see that finding a set of orthogonal measurement vectors to maximize 
the probability of detection is equivalent to the maximization problem of (p2|)-(p3|). 



Necessary and sufficient conditions for an optimum measurement maximizing (84) have been 



derived [p9| , 4C, 38|. However, except in some particular cases [38, 41, 42], obtaining a closed- form 
analytical expression for the optimal measurement directly from these conditions is a difficult and 



unsolved problem. Iterative algorithms for maximizing (84) for arbitrary vector sets are given in 

MM- 

We conclude that in general there is no known analytical expression for the GULSF. In practice. 



the GULSF may be obtained using the iterative algorithms of |^, |33]- A more detailed discussion 
on the GULSF can be found in |37]. 



9 Distance Properties of GU Frames 

So far we have mainly been concerned with structural properties of GU frames. In this section we 
study the Euclidean distance properties of GU frames. 

Suppose we are given a GU frame {(j)i = Ui(p,l < i < n} generated by the group Q with 
= 1. We would like to characterize the distance profile a{i,j) = \ \(pi — (pjW^ for all 

Since the vectors (pi are geometrically uniform, {a{i,j),l < j < n} is just a permutation of 
{d{i) = 110 — < i < n}. Furthermore, 



d{{) = {(/>, 0) + {Ui<P, Uict)) - 2^{{UicP, <P)) = 2(1 - 5R(ai)), 



(85) 
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where {oj, 1 < i < n} are the elements of the first row of the Gram matrix corresponding to the 
frame 

In applications it may be desirable to construct a GU frame such that d{i) > for 2 < i < n. 
Since d{i) = \ — Ui)(p\\'^, a sufficient condition is that 

det(/ -Ui)j^O, 2<i<n, (86) 

which is satisfied if and only if none of the matrices Ui has an eigenvalue equal to 1. Note that if 
( |86|) is satisfied, then d{i) > 0, 2 < i < n regardless of the generating vector (j) ^ 0. Groups with 
unitary representations satisfying (|8^) are known as fixed-point free groups, and have been studied 
extensively in the literature (see e.g., [H]). Thus if Q is a representation of a fixed-point free group, 
then we have that d{i) > 0,2 < i < n. 

Fixed-point free groups have recently been studied in the context of unitary space-time codes 
[^ ]. In particular, it was shown in (3^ that an abelian group of matrices {Ui} satisfies (|86| ) if and 
only if it is cyclic, i.e., Ui = W with [/" = I, and where U can be parameterized as 

U = diag(e^'2™i/'^, . . . , (87) 

where is relatively prime to n for all k. An optimization over the Uk can be performed to obtain 
distance profiles with certain prescribed properties (there are </>(n) positive integers less than n that 
are relatively prime to n, where (p{n) denotes the Euler totient function of n). 
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10 Conclusion 



In this paper we introduced the concept of GU and CGU frames and discussed some of their key 
properties. These frames may be viewed as generalizations of WH frames and filter bank frames to 
the group case. A fundamental characteristic of these frames is that they posses strong symmetry 
properties that may be desirable in a variety of applications. In particular, like WH frames and 
wavelet frames, GU frames are generated by a single generating vector. Furthermore, the canonical 
and dual frame vectors associated with a GU frame are themselves GU and are therefore also 
generated by a single generating vector which can be computed very efficiently using a Fourier 
transform matrix defined over an appropriate group. 

We also showed that GU frame vectors posses interesting symmetry properties when one or 
more frame elements are removed. This property of GU frames may be of importance in multiple 
description source-coding where it is often desirable that the quality of the reconstruction should 
not depend on the particular element lost (removed). 

Although in this paper we have focused on the case in which the underlying group is a finite 
abelian group, many of the results can be extended to the more general case of infinite-dimensional 
and non-abelian groups. An interesting direction for further research is to characterize these more 
general cases of GU frames using possibly continuous-time Fourier transforms defined over non- 
abelian groups (see e.g., [^]). 
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